Dynamics and Manipulation of Matter- Wave Solitons in Optical Superlattices 
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We analyze the existence and stability of bright, dark, and gap matter- wave solitons in optical 
superlattices. Then, using these properties, we show that (time-dependent) "dynamical superlat- 
tices" can be used to controllably place, guide, and manipulate these solitons. In particular, we 
use numerical experiments to displace solitons by turning on a secondary lattice structure, transfer 
solitons from one location to another by shifting one superlattice substructure relative to the other, 
and implement solitonic "path- following" , in which a matter wave follows the time-dependent lattice 
substructure into oscillatory motion. 
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INTRODUCTION 



After the first experimental realization of Bose- 
Einstein condensates (BECs) in dilute alkali metal va- 
pors Q, their study has experienced enormous exper- 
imental and theoretical advancements Their po- 
tential applications — ranging from matter-wave optics 
to precision measurements and quantum information 
processing — are widely held to be very promising. 

External electromagnetic fields or laser beams are used 
to produce, trap, and manipulate BECs. Additionally, 
using highly anisotropic traps, it is possible to produce 
quasi one-dimensional (ID) BECs (see, e.g., 0). In the 
formation of the latter (lying, e.g., along the x-direction), 
atoms are trapped using a confining magnetic or opti- 
cal potential V{x). In early experiments, only harmonic 
potentials were employed, but a wide variety of poten- 
tials can now be implemented experimentally. Among 
the most frequently studied, both experimentally and 
theoretically (see, e.g., |^ and references therein), are 
periodic optical lattice potentials created using counter- 
propagating laser beams Such potentials have been 
used to study Josephson effects 'g^, squeezed states 0, 
Landau-Zener tunneling and Bloch oscillations |^, the 
classical 9] and quantum 10] superfluid-Mott insulator 
transitions, and so on. Additionally, with each lattice site 
occupied by one alkali atom in its ground state, BECs in 
optical lattices show promise as registers for quantum 
computers llj. Optical lattice potentials are, therefore, 
of particular interest from the perspective of both fun- 
damental quantum physics and its connection to appli- 
cations. 

An important generalization of this setting was re- 
cently realized experimentally when ^^Rb atoms were 
loaded into an optical "superlattice" by the sequential 
creation of two lattice structures 12]. Stationary super- 



lattices can be described mathematically in the form 

V{x) = Vi cos{f<iix) + V2 cos{f<i2x) , (1) 

where a^i and i<i2 > 1^1 are, respectively, the primary and 
secondary lattice wave numbers, and Vi and V2 are the 
associated sublattice amplitudes. The great flexibility 
of superlattice potentials arises from the fact that the 
above parameters can be tuned experimentally, provid- 
ing precise control over the shape and time- variation of 
the external potential. Nevertheless, despite the afore- 
mentioned experiments, there have thus far been very 
few theoretical studies of BECs in superlattice potentials; 
these include work on dark and gap 0| solitons, the 
Mott-Peierls transition non- mean- field effects fl^ . 
and spatially extended solutions 0, . 

The aim of this work is to show that optical superlat- 
tice potentials may be used not only to sustain solitary 
matter-waves but also to manipulate them at will. As 
we illustrate below, the addition of the secondary lat- 
tice makes optical superlattices considerably more flex- 
ible than regular optical lattices. Using effectively ID 
settings, our model is the ID Gross-Pitaevskii (CP) equa- 
tion PJI in the following dimensionless form , 
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■ ^ g\u\^u^V{x)u, 



(2) 



where u is the mean-field BEC wavefunction, the non- 
linearity coefficient ^ = ±1 accounts for repulsive and 
attractive interatomic interactions, respectively, and the 
potential V{x) is given by Eq. In this paper, we 

study the kinematics, stability, and dynamics of bright 
solitons (for g = —1) as well as dark and gap solitons 
(for ^ = +1). We subsequently utilize "dynamical super- 
lattices" (in which specific lattice parameters are time- 
dependent) to show that superlattice potentials can be 
used to controllably guide, deposit, and manipulate soli- 
tons. Because of this flexibility, matter- wave solitons 
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FIG. 1: The top panel shows the branch of stable solitons 
(centered at xq ^ 2.5), with the normalized number of parti- 
cles N — as a function of the frequency A. The 
branch terminates exactly at the edge of the first band of 
excitations (the shaded rectangle) of the underlying linear 
spectrum. The bottom panel shows (for A = — 1) the profile 
of the solution (thick solid curve), together with the super- 
lattice potential (dashed curve), the primary lattice potential 
(dotted curve), and the effective potential felt by the soliton 
(soUd curve). 

(employed as information carriers) loaded into superlat- 
tice potentials may prove useful for quantum computing 
applications. 

BRIGHT SOLITONS 

For attractive interactions {g = —1), and in the ab- 
sence of any potential \y{x) = 0], Eq. Q possesses an 
exact stationary bright soliton solution of the form, 

u{x^ t) = T] sech [r]{x — xq)] exp(— zAt), (3) 

where rj is the amplitude (and inverse width), A = 
—r]^/2 = /i is the frequency (i.e., the chemical potential), 
and Xq is the location of the soliton center. Traveling 
solitons with a constant velocity can also be generated 
by applying a Galilean boost to the solution in Eq. Q. 
Bright solitons have been realized in experiments 20], 
and it is also feasible to generate them in optical lattices 
and superlattices. 

In the presence of the potential, Eq. is a perturbed 
Hamiltonian system with perturbation energy 

/oo 
V{x)\u\'^dx. (4) 
-CO 

The reduced Hamiltonian obtained by inserting the 
solutions into t he p erturbation energy is readily 
evaluated to be HHilQ 

H{xo)= V . uT^^^/o \ cos{i^jXo) . (5) 



According to Refs. 0, [23, E^, stationary states of the 
perturbed system can be obtained by demanding that 
dH{xo)/dxo = 0. That is, the selected center positions 
xo are critical points of the reduced Hamiltonian. Fol- 
lowing this perturbative treatment, one can also derive 
the equation of motion describing the center of the soli- 
ton as if it were a particle in an effective potential |2^ : 

d'^xp _ 1 dH{xo) 1 dVef^jxp) 
dt'^ T] dxp 2r] dxp 

Furthermore, as shown rigorously in 23], the small, 
nonzero linear stability eigenvalues (of the linearization 
around the solitary wave) in the presence of the pertur- 
bation are given by = —irid^H{xp)/dxp. Hence, the 
soliton is stable (unstable) at the minima (maxima) of 
the effective potential. If initialized at a maximum, the 
soliton can be seen to split (a) symmetrically under time- 
evolution of Eq. (0) if the two minima next to the max- 
imum have the same height and (b) asymmetrically if 
they have different heights. 

The top panel of Fig. ^ shows the parametric con- 
tinuation of the (always) stable bright soliton solution 
centered at the minimum of the effective potential (see 
the bottom panel). We have also obtained solutions cen- 
tered at the maxima of the effective potential, but they 
are always dynamically unstable (and, hence, are not 
shown here). We have tried to identify solutions cen- 
tered at intermediate points, as perhaps suggested by 
the metastable maxima and minima of the regular po- 
tential V{x) (see the bottom panel of Fig. The fact 
that this has not been possible — the effective potential 
lacks the metastable critical points of V{x) (again, see 
the bottom panel of Fig. — confirms that the effective 
potential^ rather than the actual one, is governing the 
steady states. Unless otherwise stated, we used the pa- 
rameter values Vi = V2 = 0.5 and tZ2 = 3/^i = 3.75 
in our numerical simulations. The ratio 1^2/1^1 = 3 was 
chosen to correspond to the superlattice reported experi- 
mentally in Ref. [l^ . We consider Vi = V2 for simplicity, 
but similar results can be obtained with other sublattice 
amplitudes. 

The linearization of Eq. Q describes the case of a non- 
interacting condensate (^ = 0). Its solution is given by a 
superposition of Bloch waves and its spectrum consists of 
bands of eigenvalues (frequencies) A = A^(/c), where the 
index n denotes the band index and the quasimomentum 
/c is a real wavenumber of bounded Bloch waves [2^ . Dif- 
ferent bands are separated fw^'gaps" in which lm{k) 7^ 0. 
As discussed in, e.g., Refs. [jjlll^? ^ superlattice poten- 
tial yields a more complicated linear band structure than 
a regular lattice, as it includes smaller gaps (called "mini- 
gaps" in [13, 13) addition to regular gaps. With the 
above parameter choices, the underlying linear spectrum 
for the superlattice potential has its first band in the in- 
terval [—0.1643, —0.1178] (see the shaded area in the top 
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FIG. 2: Unstable dark soliton centered at a potential min- 
imum (xo ~ 2.5) with A = 1. The top left panel shows 
the stationary spatial profile u{x) (soUd) and the super lat- 
tice potential (dashed). The spectral plane (top right) of 
the Unearization eigenvalues A = Ar + iXi reveals a so-called 
"loxodromic quartet" [i.e., all four eigenvalues have the same 
(|Re(A)|, |Im(A)|)] that causes the instability. The bottom 
panel, showing the magnitude of the real part of the eigen- 
value quartet vs. the frequency A, indicates that this soliton 
branch becomes stable for A < 0.55. 
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FIG. 3: A branch of gap solitons in the finite gap between the 
first and second bands. The top panel shows the normalized 
number of particles N in each soliton and the lower two panels 
show a typical example of the solution (around the middle 
of the branch) for A = 0.15. The stability of the solitons 
depends on the value of the chemical potential (and may also 
be affected for finite domains by the size of the computational 
domain) . 



panel of Fig.QJ, its second band in [0.39, 1.5244] (with a 
minigap in [0.7748,0.8171]), and so on. 



DARK AND GAP SOLITONS 



For the superlattice potential (QJ, Eq. Q becomes 



(Pxq 1 1 

= sin(/^iXo) + -f^2V2 sin(A^2^o) • (10) 



We now turn to the study of repulsive BECs with g = 1 
and first consider dark soliton solutions of Eq. Dark 
solitons have been realized experimentally 25] in BECs 
confined in harmonic traps and have been studied theo- 
retically in lattices [2^ [23 and superlattices p^] , where 
they can also be created. 

We seek solutions of Eq. ^ of the form u = 
0(x, t) exp(— it) (we set the chemical potential /i = 1) and 
assume that the potential, characterized by a scale i?, is 
slowly- varying on the soliton scale (which is on the order 
of the healing length). Then, following the multiple-scale 
boundary layer theory of Ref. 28] (similar results can be 
obtained using other perturbative approaches 29]), we 
use the ansatz 

^= [^o + ^^i(2:-xo,t)]exp[zi9 + z^(x-xo)]+0(£:^) (7) 

for the interval |x — xo| < where 

(j)o = i{vo — v) -\- /ctanh[/c(x — xq)] , (8) 

with vo = xo and 0{t) = (1/2) [0{xo - R,t) ^ 0{xo + R, t)] 
{v is defined analogously). Letting the cut-off R ^ oo 
and considering solutions satisfying v = dO/dx = yields 
the equation of motion for the soliton center xq, 

d^xp _ 1 dV{xo) dVefi . . 

df^ ~ 2 dxo ~ dxo ' ^ ^ 



Figure 121 shows a dark soliton centered at a minimum 
of the superlattice. The depicted solution is unstable 
due to a loxodromic quartet of eigenvalues but can be- 
come stable for a frequency A < 0.55 (see the bottom 
panel in Fig. [2|). A solution centered at a maximum of 
the potential is unstable due to the presence of a large 
real eigenvalue pair (e.g., for A = 1, the eigenvalues have 
magnitude 0.3833). We have also attempted to identify 
solutions centered between consecutive absolute maxima 
and minima, i.e., near metastable minima and maxima 
of the effective potential (e.g., for x ^ 0.9 and x ~ 1.6), 
but were unable to find any. This apparent violation of 
the effective energy landscape suggested by the bound- 
ary layer theory poses an interesting question for future 
studies. 

In the same context of repulsive condensates (again 
with ^ = 1), we have also obtained gap soliton solutions 
of Eq. ((21). Such gap solitons are spatially localized non- 
linear modes that occur in the band gaps of the linearized 
spectrum and, in fact, have the form of bright solitons but 
in a repulsive medium. Figure [2l shows a stable branch of 
gap solitons that exists throughout the entire finite gap 
between the first and second bands. Gap solitons have 
recently been obtained experimentally in regular optical 
lattices [s^ , and thus they can also be straightforwardly 
produced experimentally in the superlattice setting. 
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Bright soliton transfer 



FIG. 4: (Color Online). Bright soliton evolution for the 
V2 cos{Hi2x) potential. The soliton is initially centered at 
the minimum near 0.8, which becomes a met ast able mini- 
mum of the super lattice potential V{x) and disappears alto- 
gether as a minimum for the effective potential Veff{x) once 
the Vi cos{k,ix) potential is turned on. The abrupt switching- 
on of the second sublattice at t = 20 is represented function- 
ally by Vi{t) = (1/2)[1 +tanh(5(t - 20))]. The soliton can 
no longer stay in its original location, so it goes to the clos- 
est minimum of the effective potential (solid red/light curve), 
about which it oscillates. 



DYNAMICAL SUPERLATTICES 

We now turn to using the superlattice as a means to 
guide, displace, and (more generally) manipulate matter- 
wave solitons at will in the potential landscapes discussed 
above. In recent experiments 31], the center of a regu- 
lar optical lattice was "shaken" (translated periodically) 
to examine a period-doubling instability in BECs. Such 
dynamical manipulations of superlattice potentials are 
similarly achievable in the laboratory. 

Given a stable bright soliton located at a minimum of a 
regular optical lattice, we can displace it from its location 
by turning on the second sublattice abruptly (nonadia- 
batically), as indicated by the space-time plot in Fig. ^ 
As the soliton is no longer located at a minimum of the 
effective potential, it cannot stay in its original location, 
so it moves to the closest minimum and oscillates in that 
well. 

As another example, in Fig. O we show how to use a 
sublattice of the superlattice as a means of transferring 
bright and gap solitons at will. A soliton is placed in 
the potential Vi cos{Hii{x — Xc)) + V2 cos(/t:2(^)), the first 
sublattice of which is subsequently displaced according 
to the following equation. 



tanh 



t-to 



(11) 



where Xi and Xf denote, respectively, the initial and final 
"center" positions and the parameter r determines the 
speed of the displacement 113 • We considered both 
adiabatic cases (large r; left panels), where the soliton 
transfer can be successful, and nonadiabatic ones (small 
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FIG. 5: (Color Online). Controllable transfer of bright (top 
quartet) and gap (bottom quartet) solitons by manipulating 
one sublattice of the superlattice potential. For each quartet, 
the top panels show space-time plots (with colors indicat- 
ing the value of |i^|^) and the bottom panels show the time- 
evolution of the soliton center (solid) and superlattice "cen- 
ter" Xc (dashed). The left panels show the result of adiabatic 
potential displacement (r = 5) and the right panels show 
the result of nonadiabatic potential displacement (r = 0.5). 
Clearly, adiabaticity plays a major role in the success of the 
soliton transfer. 



r; right panels), where such manipulations will typically 
fail to guide the soliton. Note that similar results have 
also been obtained both for bright and dark soli- 
tons but in regular optical lattices. 

We also performed a more "demanding" experiment, 
in which the bright and gap solitons were not merely 
deposited at a new location but were instead "instructed" 
to follow the time-dependent sublattice into oscillatory 
motion. This numerical experiment, shown in Fig. |S1 
corresponds to a lattice displaced according to 



Xf , (t-to 

Xc = Xi^ ^ sm 

2 V Ti 



1 + tanh 



(12) 



We turn on the sinusoidal potential abruptly (with r = 
0.5), but vary its period (using both large and small 
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Bright soliton path- following 
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FIG. 6: (Color Online). As in Fig.|^ but for controllable path- 
following of bright (top quartet) and gap (bottom quartet) 
solitons. The evolution of the soliton center is governed by 
Eq. (fT^ for n = 10 (left panels) and n = 5 (right panels) 
[t = 0.5 in both cases]. Clearly, it is necessary that the lattice 
oscillations be sufficiently adiabatic for solitons to follow a 
dynamical superlattice's path successfully. 



Ti). We observe that while the solitons lose power (i.e., 
f dx\u\'^) in traversing the rough terrain of the immobile 
lattice (in this example, the bright soliton with the larger 
power emits more radiation) , they can still follow the os- 
cillation, provided the motion is, again, sufficiently adi- 
abatic. For small ri (the nonadiabatic case), the waves 
completely disintegrate into small-amplitude radiation. 



CONCLUSIONS 

In this work, we analyzed the properties of bright, 
dark, and gap matter-wave solitons in the presence of 
superlattice potentials. We focused, in particular, on 
showing (for typical potential parameter values) how the 
dynamical modification of the (experimentally tunable) 
properties of the superlattice potential can be used to 
guide, transfer, and deposit these coherent structures 



across their corresponding energy landscape. We be- 
lieve that such controllable soliton manipulation bene- 
fits greatly from the enhanced flexibility arising from the 
extra length scale and tunable parameters of the super- 
lattice structure as compared to regular optical lattices. 
This, in turn, may pave the way towards the macroscopic 
manipulation of solitonic "bits" of information in BECs 
in a manner that bears similarities to (but also enables 
extensions of) the setting of nonlinear optics. 
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